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The LHCb Collaboration has recently observed evidence of CP vio- 
lation in neutral D meson decays. CP violation in the charm sector is 
generically expected to be very small in the Standard Model, but can be 
enhanced in many models of new physics. In this document we will present 
the results of a search for time-integrated CP violation in h~h~^ with 

{h = K,7r) decays, performed with around 0.6 fb^^ of data collected by 
LHCb in 2011. The difference in CP asymmetry between K~K^ 
and — >■ tt^tt"*", ^Aqp = Acp{K~K^) — Acpi'K'T^^) is measured to 
be AAcp = [-0.82 ± 0.21(stot.) ± 0.11(syst.)]% . This differs from the 
hypothesis of CP conservation by 3.5 sigma. 
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1 Introduction 



In the past few years, thanks to the observation of mixing [T|[2] the interest for 
charm physics has increased significantly. Mixing is now well established [sj at a level 
which is consistent with expectations |4] . Recently LHCb reports the first evidence 
of CP violation in singly Cabibbo-suppressed decays of mesons to two-body final 
states. This result has definitely heightened the theoretical interest in charm physics. 
Before LHCb measurement, CP asymmetries in these decays were expected to be very 
small in the SM (5]-[8], with naive predictions of up to O(10~'^). In this scenario, the 
LHCb result, characterised by a central value of O(10~^), has been unexpected. In 
this sector, it is very difficult to achieve precise theoretical predictions of CP violation. 
This is due to the fact that the charm quark is too heavy for chiral perturbation and 
too light for the heavy-quark effective theory to be applied reliably. This leads to 
no conclusion if the observed effect is a clear sign of New Physics [9 12 . In order 



to clarify this scenario, it is necessary to carry out measurements in other charm 
decays. In these proceedings, LHCb results for the measurement of the difference in 
integrated CP asymmetries between — t- K~ and — )■ tt^tt"'", performed with 
0.6 fb~^ of data collected in 2011, are presented. 



2 The LHCb detector 
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The LHCb experiment [13j is designed to exploit the huge hh cross section 
pp collisions LHC energies. However thanks to the large charm cross section of 
(6.10±0.93) mb in 7 TeV proton-proton collisions fls], LHCb can also make excellent 
measurement in the charm sector. The LHCb detector is a single arm spectrometer 
in the forward direction. It is composed of a vertex detector around the interaction 
region, a set of tracking stations in front of and behind a dipole magnet that provides 
a field integral of 4 Tm, two Ring-Imaging Cherenkov (RICH) detectors, electromag- 
netic and hadronic calorimeters complemented with pre-shower and scintillating pad 
detectors, and a set of muon chambers. The two RICH detectors are of particular 
importance for this analysis, as they provide the particle identification (PID) infor- 
mation needed to disentangle the various D — )• h^h'~ final states. They are able to 
separate efficiently vr, K and in a momentum range from 2 GeV/ c up to and beyond 
100 GeV/ c. RICH-1 is installed in front of the magnet and uses areogel and C4F10 
as radiators, while RICH-2 is installed behind the magnet and employs CF4. 
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3 Formalism 



The time-dependent CP asymmetry Acp{f'., t) for decays to a CP eigenstate / 
(with f = f) is defined as 

^ rW)^/)-r(g°W^/)^ 

where F is the decay rate for the process indicated. In general Acp{f] t) depends 
on /. For / = K~K^ and / = tt^tt^, Acp{f'i t) can be expressed in terms of 
two contributions: a direct component associated with CP violation in the decay 
amplitudes, and an indirect component associated with CP violation in the mixing 
or in the interference between mixing and decay. The asymmetry Acp{f'i t) may be 



written to first order as 16, 17 



Acp{f; t) = a^^(/) + ^a-^, (2) 

where cicp(/) is the direct CP asymmetry, r is the lifetime, and a'^p is the indirect 
CP asymmetry. To a good approximation this latter quantity is universal |7,18 . The 
time-integrated asymmetry measured by an experiment, Acp{f), depends upon the 
time-acceptance of that experiment. It can be written as 



AcpU) = a'^Jpif) + ^a-;^, (3) 

where (t) is the average decay time in the reconstructed sample. Denoting by A the 
differences between quantities for — )■ and — t- 7r^7r+ it is then possible 

to write 

A^cp = Acp{K-K+) - Acp(7r-7r+) (4) 

r^dir / 7>— ^dir/ , ^(^) „md 

In the limit that A(t) vanishes, AAcp is equal to the difference in the direct CP 
asymmetry between the two decays. However, if the time-acceptance is different for 
the K~K^ and Tr'Tr"*" final states, a contribution from indirect CP violation remains. 

What LHCb measures is the raw asymmetry for tagged decays to a final state 
/ is given by Araw(/)- It is defined as 

where N{X) refers to the number of reconstructed events of decay X after background 
subtraction. 
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To first order the raw asymmetries may be written as a sum of four components, 
due to physics and detector effects: 



Aaw(/) = Acp{f) + Auif) + AD(7r+) + Ap{D*+). 



(6) 



Here, ^d(/) is the asymmetry in selecting the decay into the final state /, ^0(71+) 
is the asymmetry in selecting the slow pion from the D*^ decay chain, and Ap{D*^) is 
the production asymmetry for D*~^ mesons. The asymmetries and Ap are defined 
in the same fashion as Araw The first-order expansion is valid since the individual 
asymmetries are small. 

For a two-body decay of a spin-0 particle to a self-conjugate final state there can 
be no detection asymmetry, i.e. Ap,{K~ K^) = Aj){'n-^7i^) = 0. Moreover, Aj){7i'^) 
and Ap{D*~^) arc independent of / and thus in the first-order expansion of equation 
5 those terms cancel in the difference Ay.aw{K^ K^) — ylraw(7r~7r+), resulting in 



To minimise second-order effects that are related to the slightly different kinematic 
properties of the two decay modes and that do not cancel in AAcp, the analysis is 
performed in bins of the relevant kinematic variables, as discussed later. 

4 Event selection 

Selections are applied to provide samples of D*^ — > -D°vr+ candidates, with D° — )■ 
K~ or 7r~7r"'". Events are required to pass both hardware and software trigger 
levels. A loose D° selection is apphed in the final state of the software trigger, and 
in the offline analysis only candidates that are accepted by this trigger algorithm are 
considered. Both the trigger and offline selections impose a variety of requirements 
on kinematics and decay time to isolate the decays of interest, including requirements 
on the track fit quality, on the and D*^ vertex fit quality, on the transverse 
momentum (px > 2 GeV/c) and decay time {ct > 100 |^m) of the D° candidate, on 
the angle between the D° momentum in the lab frame and its daughter momenta in 
the D'^ rest frame (|cos^| < 0.9), that the D° trajectory points back to a primary 
vertex, and that the daughter tracks do not. In addition, the offline analysis 
exploits the capabilities of the RICH system to distinguish between pions and kaons 
when reconstructing the D'^ meson, with no tracks appearing as both pion and kaon 
candidates. 

A fiducial region is implemented by imposing the requirement that the slow pion 
lies within the central part of the detector acceptance. This is necessary because the 
magnetic field bends pions of one charge to the left and those of the other charge to 
the right. For soft tracks at large angles in the xz plane this implies that one charge 



AAcP^A,^^{K-K+) - Aaw(7r-7r+). 



(7) 
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Figure 1: Fits to the (a) m{K~ K'^) and (b) m(7r^7r+) spectra of D*^ candidates 
passing the selection and satisfying < 5m < 15 MeV/c^. The dashed hne corre- 
sponds to the background component in the fit, and the vertical lines indicate the 
signal window of 1844-1884 MeV/c^. 

is much more likely to remain within the 300 mrad horizontal detector acceptance, 
thus making ^^(Tr^) large. Although this asymmetry is formally independent of 
the decay mode, it breaks the assumption that the raw asymmetries are small 
and it carries a risk of second-order systematic effects if the ratio of efficiencies of 

— )■ K^K^ and — t- 7r^7r+ varies in the affected region. The fiducial requirements 
therefore exclude edge regions in the slow pion {px,p) plane. Similarly, a small region 
of phase space in which one charge of slow pion is more likely to be swept into 
the beam-pipe region in the downstream tracking stations, and hence has reduced 
efficiency, is also excluded. After the implementation of the fiducial requirements 
about 70% of the events are retained. 

The invariant mass spectra of selected K~K^ and ti~tt^ pairs are shown in Fig. [!} 
The half-width at half-maximum of the signal line-shape is 8.6 MeV/c^ for K~K^ 
and 11.2 MeV/c^ for h'tt^ , where the difference is due to the kinematics of the decays 
and has no relevance for the subsequent analysis. The mass difference {5m) spectra 
of selected candidates, where 5m = m{h^h'^'K^) — m{h^h^) — m(7r+) for h = KjTt, 
are shown in Fig. |2] 

Candidates are required to lie inside a wide 5m window of 0-15 MeV/c^, and in 
Fig. |2] and for all subsequent results candidates are in addition required to lie in a 
mass signal window of 1844-1884 MeV/c^. The D*~^ signal yields are approximately 
1.44 X 10^ in the K~K^ sample, and 0.38 x 10^ in the tt~h^ sample. Charm from h- 
hadron decays is strongly suppressed by the requirement that the originate from a 
primary vertex, and accounts for only 3% of the total yield. Of the events that contain 
at least one D*^ candidate, 12% contain more than one candidate; this is expected 
due to background soft pions from the primary vertex and all candidates are accepted. 
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Figure 2: Fits to the 5m spectra, where the D*^ is reconstructed in the final states 
(a) K^K^ and (b) tt^tt^, with mass lying in the window of 1844-1884 MeV/c^. The 
dashed line corresponds to the background component in the fit. 

The background-subtracted average decay time of candidates passing the selection 
is measured for each final state, and the fractional difference A{t)/T is obtained. 
Systematic uncertainties on this quantity are assigned for the uncertainty on the 
world average lifetime r (0.04%), charm from 6-hadron decays (0.18%), and the 
background-subtraction procedure (0.04%). Combining the systematic uncertainties 
in quadrature, we obtain A(t)/r = [9.83 ± 0.22(stat.) ± 0.19(syst.)] %. The 7r-7r+ 
and K~K^ average decay time is (t) = (0.8539 ± 0.0005) ps, where the error is 
statistical only. 

5 Fit procedure 

Fits are performed on the samples in order to determine Aj.i^„{K~ K'^) and y4raw(7r^7r+). 
The production and detection asymmetries can vary with and pseudo-rapidity rj, 
and so can the detection efficiency of the two different decays, in particular through 
the effects of the particle identification requirements. The analysis is performed in 
54 kinematic bins defined by the pt and rj of the D*~^ candidates, the momentum 
of the slow pion, and the sign of px of the slow pion at the D*^ vertex. The events 
are further partitioned in two ways. First, the data are divided between the two 
dipole magnet polarities. Second, the first 60% of data are processed separately from 
the remainder, with the division aligned with a break in data taking due to an LHC 
technical stop. In total, 216 statistically independent measurements are considered 
for each decay mode. 

In each bin, one-dimensional unbinned maximum likelihood fits to the Sm spectra 
are performed. The signal is described as the sum of two Gaussian functions with 
a common mean /i but different widths cXj, convolved with a function B{6m;s) = 
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Q{Sm) 6m'^ taking account of the asymmetric shape of the measured 6m distribution. 
Here, s ~ —0.975 is a shape parameter fixed to the value determined from the global 
fits shown in Fig. |2| B is the Heaviside step function, and the convolution runs 
over 6m. The background is described by an empirical function of the form 1 — 
g-{(5m-(5m,o)/a^ where Smo and a are free parameters describing the threshold and 
shape of the function, respectively. The D*~^ and D*~ samples in a given bin are 
fitted simultaneously and share all shape parameters, except for a charge-dependent 
offset in the central value fi and an overall scale factor in the mass resolution. The 
raw asymmetry in the signal yields is extracted directly from this simultaneous fit. 
No fit parameters are shared between the 216 subsamples of data, nor between the 
K^K^ and vr^vr"*' final states. 

The fits do not distinguish between the signal and backgrounds that peak in 6m. 
Such backgrounds can arise from D*~^ decays in which the correct slow pion is found 
but the is partially mis-reconstructed. These backgrounds are suppressed by 
the use of tight particle identification requirements and a narrow mass window. 
From studies of the D° mass sidebands (1820-1840 and 1890-1910 MeV/c^), this 
contamination is found to be approximately 1% of the signal yield and to have small 
raw asymmetry (consistent with zero asymmetry difference between the K~K^ and 
TT~T[^ final states). Its effect on the measurement is estimated in an ensemble of 
simulated experiments and found to be negligible; a systematic uncertainty is assigned 
below based on the statistical precision of the estimate. 

6 Result 

A value of AAcp is determined in each measurement bin as the difference between 
Aj-3_y^{K~ K^) and y4i.aw(7r^7r+). Testing these 216 measurements for mutual consis- 
tency, we obtain x^/ndf = 211/215 (x^ probability of 56%). A weighted average 
is performed to yield the result AAcp = (—0.82 ± 0.21)%, where the uncertainty is 
statistical only. 

Numerous robustness checks are made. The value of AAcp is studied as a function 
of the time at which the data were taken (Fig. [s]) and found to be consistent with 
a constant value (x^ probability of 57%). The measurement is repeated with pro- 
gressively more restrictive RICH particle identification requirements, finding values 
of (-0.88 ± 0.26)% and (-1.03 ± 0.31)%; both of these values are consistent with 
the baseline result when correlations are taken into account. Other checks include 
applying electron and muon vetoes to the slow pion and to the daughters, use of 
different kinematic binnings, validation of the size of the statistical uncertainties with 
Monte Carlo pseudo-experiments, tightening of kinematic requirements, testing for 
variation of the result with the multiplicity of tracks and of primary vertices in the 
event, use of other signal and background parameterizations in the fit, and imposing 
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Figure 3: Time-dependence of the measurement. The data are divided into 19 disjoint, 
contiguous, time-ordered blocks and the value of AAcp measured in each block. The 
horizontal red dashed line shows the result for the combined sample. The vertical 
dashed line indicates the technical stop. 



a full set of common shape parameters between D*~^ and D*~ candidates. Poten- 
tial biases due to the inclusive hardware trigger selection are investigated with the 
subsample of data in which one of the signal final-state tracks is directly responsible 
for the hardware trigger decision. In all cases good stability is observed. For sev- 
eral of these checks, a reduced number of kinematic bins are used for simplicity. No 
systematic dependence of AAcp is observed with respect to the kinematic variables. 



7 Systematic uncertainties 

Systematic uncertainties are assigned by: loosening the fiducial requirement on the 
slow pion; assessing the effect of potential peaking backgrounds in Monte Carlo 
pseudo-experiments; repeating the analysis with the asymmetry extracted through 
sideband subtraction in 6m instead of a fit; removing all candidates but one (chosen 
at random) in events with multiple candidates; and comparing with the result ob- 
tained without kinematic binning. In each case the full value of the change in result 
is taken as the systematic uncertainty. These uncertainties are listed in Table [1} The 
sum in quadrature is 0.11%. Combining statistical and systematic uncertainties in 
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Table 1: Summary of absolute systematic uncertainties for AAcp- 



Source Uncertainty 

Fiducial requirement 0.01% 

Peaking background asymmetry 0.04% 

Fit procedure 0.08% 

Multiple candidates 0.06% 

Kinematic binning 0.02% 

Total 0.11% 



quadrature, this result is consistent at the la level with the current HFAG world 
average [s]. 

8 Conclusion 

The time-integrated difference in CP asymmetry between K^K^ and — ^ 

TT^TT"^ decays has been measured to be 

AAcp = [-0.82 ± 0.21(stat.) ± 0.11(syst.)] % 

with 0.62 fb^^ of 2011 data. Given the dependence of AAcp on the direct and 
indirect CP asymmetries, shown in Eq. Q, and the measured value A(t)/r = 
[9.83 ± 0.22(stat.) ± 0.19(syst.)] %, the contribution from indirect CP violation is sup- 
pressed and AAcp is primarily sensitive to direct CP violation. Dividing the central 
value by the sum in quadrature of the statistical and systematic uncertainties, the 
significance of the measured deviation from zero is 3.5a. This is the first evidence 
for CP violation in the charm sector. To establish whether this result is consistent 
with the SM will require the analysis of more data, as well as improved theoretical 
understanding. 
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